Dynamic Model of Fiber Bundles 
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A realistic continuous-time dynamics for fiber bundles is introduced and studied both analytically 
and numerically. The equation of motion reproduces known stationary-state results in the deter- 
ministic limit while the system under non-vanishing stress always breaks down in the presence of 
noise. Revealed in particular is the characteristic time evolution that the system tends to resist the 
stress for considerable time, followed by sudden complete rupture. The critical stress beyond which 
the complete rupture emerges is also obtained. 
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During past decades fiber bundle models have received 
considerable attention and been studied extensively 0. 
Originally introduced to explain ruptures in heteroge- 
neous material under tension fiber bundle models 
have been applied to cracks and fractures^arthquakes, 
and other related breakdown phenomena U 0, 0| . The 
main issue in the studies of fiber bundles is the life time 
of a bundle under given applied stress and accordingly, 
how the system fails as the applied stress is gradually 
increased jfj|. Recent works have been extended to the 
study of thermodynamic transitions Q, damping effects 
of viscous fibers @, and thermally activated failures 
The breaking dynamics in the fiber bundle model is es- 
sentially determined by the condition whether the thresh- 
old of a particular fiber is exceeded by the applied stress 
and the equation of motion is expressed as recursion rela- 
tions. It is thus deterministic, even though the thresholds 
of the fibers are drawn from some probability distribu- 
tions which may include various effects [T(ij |. The recur- 
sion equations can be solved by means of the graphical 
method and have been shown to be useful in classifying 
all possible threshold distributions. Here time is viewed 
as a discrete variable and it is assumed that at each time 
step all the fibers meeting the breaking condition break 
immediately. The concern is then the number of intact 
fibers in the stationary state after evolving in discrete 
time. Unlike the stationary state, however, this deter- 
ministic recursive dynamics, based on synchronous up- 
dating of the fibers at each discrete time step, does not 
provide a realistic description of the actual dynamics. In 
real systems it is obvious that time is continuous and 
there usually exists retardation. Namely, it may take 
time for a fiber to respond to a given stress and for the 
stress to get redistributed among fibers. 

This work is the first attempt toward realistic dy- 
namics of fiber bundles, incorporating those features to- 
gether with the probabilistic generalization. This kind of 
continuous-time dynamics in general yields the equations 
of motion in the form of delay-differential equations and 
was successfully applied to neural network problems . 



From this formulation, we obtain the equation of mo- 
tion for the average number of intact fibers, and find 
that the stationary solutions for simple threshold distri- 
butions are the same as those found previously in the 
recursion-equation approach. For realistic distributions, 
the equation is solved numerically to give the time evo- 
lution. It is revealed that after initial rupture the aver- 
age number of intact fibers tends to remain more or less 
constant, exhibiting plateau-like behavior as a function 
of time. In other words, even in case that the stress is 
strong enough to break all the fibers eventually, the sys- 
tem resists the stress for considerable time before com- 
plete rupture. The critical stress beyond which the com- 
plete rupture emerges is also obtained. 

We consider a bundle of N fibers pulled by an external 
force F. In the global load-sharing limit, all fibers share 
the force uniformly so that which provides stress to each 
fiber. Each fiber has its own threshold and resists the 
stress lower than the threshold, thus remaining intact. 
If the stress exceeds the threshold, however, the fiber 
becomes broken, leaving the applied stress redistributed 
among the neighboring intact fibers. For convenient rep- 
resentation, we assign a "spin" variable to each fiber in 
such a way that Sj = +1 (— 1) for the ith fiber broken 
(intact). The state of the bundle is described by the con- 
figuration of all the fibers, i.e., s = (si, S2, sjv)- The 
total number of the intact fibers is related with the av- 
erage spin s 
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and we are interested in how N_ evolves in time as well 
as its stationary value. 

The total stress on the ith fiber can then be written in 
the form 



n 



2 ' 



(2) 



where / is the stress directly due to the external force and 
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Vij represents the stress transferred from the jth fiber (in 
case that it is broken) . The breaking of the ith fiber with 
threshold hi is determined according to: 



i]i < hi => Si = -1 
rji> hi => Si = + 1, 



(3) 



which, through the use of Eq. JJJ, can be simplified as 



SiEi > 



(4) 



with the local field Ei = (rfc — hi) (1 — s)/2. This deter- 
mines the stationary configuration at zero "temperature" 
(i.e., in the deterministic limit). 

To describe the time evolution toward the stationary 
state described by Eq. (jJJ, we also take into considera- 
tion the uncertainty ("noise") present in real situations, 
which may arise from impurities and other environmental 
influences. We thus begin with the conditional probabil- 
ity that the ith fiber breaks at time t+St, given that it is 
intact at time t: 



p(si= + l > t+St\s i =-l,t;s',t-t d ) = —[I 



tanh/lE-], 
(5) 

where s' = (s[, s' 2 , ■ ■ . , s' N ) represents the configuration of 
the system at time t—t d and E[ = (r^ — /ij)(l — s')/2 is the 
local field at time t—td. Note the two time scales td and t r 
here: td denotes the time delay during which the stress is 
redistributed among fibers while the refractory period t r 
sets the relaxation time of the system. The temperature 
T = P^ 1 measures the width of the threshold region of 
the fibers or the noise level 0]: I n the deterministic 
limit (T = 0), which is our main concern, the factor 
(1 + tanh/3£"j')/2 in Eq. reduces to the step function 
9{E' i ), yielding the stationary-state condition given by 
Eq. (0J . The conditional probability that the ith fiber is 
reconnected given that it is broken at time t is obviously 
zero 01 



p(si= - 1, t+St\si= + 1, t; a', t-t d ) = 0. 



(6) 



Equations (0 and © can be combined to give 
a general expression for the conditional probability 
p(s'i, t+St\si, t; a', t— t d ), which, in the limit 5t — > oo, can 
be expressed in terms of the transition rate: 



p(s'i,t+5t\si,t;s' ,t-t d ) 

( Wi(sf,s',t-td)6t for s- = -s. 
\ 1 - Wi(si; s', t-t d )8t for = Sj, 

where the transition rate is given by 



Wi(si;s',t-t d ) 



1 
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tanh^E- 
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The behavior of the fiber bundle is then governed by 
the master equation, which describes the evolution of 



the joint probability P(a, t; s', t— t d ) that the system is 
in state s' at time t—t d and in state s at time t: 

P(s, t+St; a', t-t d ) - P(a, t\ s', t-t d ) 



E 



[p(s", t+St\s, t; a', t-td)P(a, t; a', t-t d ) 



- p(a, t+5t\a", t; a', t-t d )P(a", t; a', t-t d )]. (9) 

Thus we obtain the equation of motion in the form of a 
non-Markov master equation. Here the conditional prob- 
ability for the whole system is given by the product of 
that for each fiber 

TV 

p(a", t+St\a, t; a', t-t d ) = *[[p(s'l, t+St\a it t; a', t-t d ). 



In the limit St — > oo, Eq. takes the differential form: 
d -P(s,t;s',t-l) = -J2[wi(si;a')P(a,t;a',t-l) 

i 

- Wl (- Sl ;s')P(Ps,i;s',i-l)],(10) 



where time t has been rescaled in units of the delay 
time t d , the transition rate is given by u)j(sj;s) = 
tdWi{si\ a', t— t d ) with Wi(si] a', t— t d ) defined in Eq. JSJl, 
and FiS = (si, S2, — Sj, Sj+i, Sjv)- Then equations 
describing the time evolution of relevant physical quanti- 
ties in general assume the form of differential-difference 
equations due to the retardation in the stress redistribu- 
tion. In particular, the average spin for the kth fiber, 
mk(t) = (sfe) = J2 S s' s kP( s , t'i s ', t— 1) can be obtained 
from Eq. (|1U|) . by multiplying Sk and summing over all 
configurations: 



j t m h {t) 



- 2_,[sfcWfc(s fc ; a')P{a, t; s', t-1) 

s.s' 

- s k Wk{-s k ;s')P(F k s, t; s', t-1)} 
-2(s k w k (s k ;a')), (11) 



where ( ) denotes the average over P(a, t;a', t— 1). Eval- 
uation of the average (skWk(sk;a')), with s| = 1 noted, 
leads to 



d 

T dt mk 



1 - m k 



1 - s k 



tanh l3E' k 



(12) 



where r = t r /t d gives the relaxation time (in units of t d ). 

To proceed further, we need to specify the explicit form 
of Vij. In the simplest case of global load sharing, which 
has been mostly studied due to the analytical tractability, 
we have = r]j/N_ as well as f — F/N. Equation @ 
then leads to (1 — s)rji = 2/ and accordingly, the local 
field Ei — f — (hi /2)(1 — s). The infinite-range nature 
of such global load sharing allows one to replace EL by 
its average (E' k ) = f — (h k /2)[l — fh(t— 1)], where it has 
been noted that s' is the configuration at time t— 1, i.e., 
(S , ) = N- 1 E j (s' j ) = N- l E j m j (t-l) = m(t-l). 
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For convenience, we now rewrite Eq. i|12fl in terms of 
the average number of intact fibers at time t. Defining 
Xk = (1 — Wfc)/2 and x = N~ x x k = (1 — m)/2, we 
have, from Eq. (N_) = Nx and thus obtain from 
Eq. I|12l) the equation of motion for the average fraction 
of intact fibers: 

T j t x k{t) = -x k (t) -x k (t)tajahp\f-h h x(t-l)] (13) 

or, upon the sample average, 

T Jt* {t) = ~ ^ ts ^^f- h ^{t-W)^ ( 14 ) 

where (( )) stands for the average with respect to the 
distribution of {ftfc} and self-averaging has been assumed. 
We first examine the stationary solutions of Eq. (|13|) 

x k + x k ta,nhf3[f-h k x)}=0, (15) 

which possesses two possible solutions: x k — or 
tanh /?(/— hkx) = —1. The latter is possible only if 
/ < h k x at T = 0. In this case, the stress / is not 
large enough to break the fiber, so that we have x k = 1. 
When / > hkx, on the other hand, Xk = is the only 
solution. It is thus concluded that 

x k = 6(h k x-f) (16) 

in the deterministic limit (T = 0). The average frac- 
tion of intact fibers, given by the sample average, i.e., 
the average over the distribution of the threshold, then 
becomes 

x = J dhg(h)x = J dh6(hx-f)g(h), (17) 

where g(h) is the distribution function of the thresh- 
old {hk}- This is the self-consistent equation, the so- 
lution of which gives the stationary fraction of the sur- 
viving fibers. Note also that the stationarity condition 
x + ({x tanh/3[/— h k x])) = or 

J dhg{h)8{hx-f)[l + tanh (3(f-hx)] = (18) 

is automatically satisfied since tanh/3(/— fix) = 
0{f~hx) - e(hx-f). At T ^ 0, we have x k = 0, giv- 
ing x = as the only possible solution. 

We now briefly discuss how to solve the Eq. (|17f) . For 
a continuous distribution of the threshold, we formally 
solve Eq. (|17fl by simply performing the integration 

x = l-G(f/x), (19) 

where G(h) = j Q dh' g(h') is the cumulative distribution 
of the threshold. This equation corresponds to the fixed- 
point equation for the deterministic dynamics, except 
that in our case x is a double averaged quantity. This 
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FIG. 1: Behavior of the average fraction x of intact fibers with 
time t (in units of the delay time) for the Gaussian distribu- 
tion of ft = 1 and a = 0.2. Solid and dotted lines correspond 
to / = 0.66 and 0.68, respectively. 

equation may be solved, e.g., graphically for given g(h) 
[lj, and yields all known results, which we do not repro- 
duce here. For example, the simple bimodal distribution 
g(h) = P 8{h-h) + (1 - P )6{h-f 2 ) with < f 2 < h and 
< p < 1 leads to 

x = J dh0{hx-f)\p8{h-h) + (1 - p)S(h-f 2 )} (20) 

whose solution is given by 

( 1 for / < f 2 x 
x = < p for f 2 X < f < fix 
{ for f lX < f. 

As / is increased from zero, the first rupture occurs at 
/ — f2i yielding x = p, and the second one occurs at 
/ = Pfi if h < Pfi- For /2 > Pfi, however, there 
appears only one rupture at / = /a. In these two cases 
the critical stress / c , beyond which the system breaks 
completely, is thus given by pfi and f 2 , respectively. 

Next the time evolution of the system is explored via 
direct numerical integration of the equation of motion 
(|13fl . We consider realistic distributions of the threshold 
including the Gaussian distribution as well as the uni- 
form one. Specifically, we set r = 100 and use the time 
step A< = 0.1, mostly in a system of N = 10 3 fibers. 
These parameter values have been varied, only to give 
no appreciable difference. 

The typical behavior of the average fraction x of in- 
tact fibers is shown in Fig. ^ for the Gaussian distri- 
bution with unit mean (ft = 1) and variance a = 0.2. 
Note that for / = 0.66 a large part of the fibers remain 
intact whereas all fibers become broken for / = 0.68, in- 
dicating that the critical stress f c lies in between. Indeed 
this is the case in Fig. which displays how the critical 
stress varies with the variance of the threshold distribu- 
tion. It is of interest that after initial rupture x tends 
not to change much and exhibits a plateau as a function 
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FIG. 2: Critical stress f c versus variance a 
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FIG. 3: Behavior of the average fraction of intact fibers with 
time at T = 0.1 for / = 0.48 (solid line), 0.58 (dotted line), 
and 0.68 (dashed line). Other parameters are the same as 
those in Fig. 
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FIG. 4: Comparison of the behavior at T — (solid line) and 
0.1 (dotted line). The stress is / = 0.68 and other parameters 
are the same as those in Fig. 



stationary features found in the deterministic recursion- 
equation approach have been reproduced. In particular 
disclosed is characteristic time evolution that the system 
tends to resist the stress for considerable time, followed 
by the complete rupture occurring suddenly. This has 
interesting implications to many systems in nature such 
as biological problems, the investigation of which is left 
for further study. 

This work was supported in part by the Overhead Re- 
search Fund of SNU and in part by the Basic Research 
Program (Grant No. R01-2002-000-00285-0) of KOSEF. 



of time, which persists regardless of the details of the 
threshold distribution. This indicates that even for the 
stress strong enough to break all the fibers eventually, 
the system resists the stress for considerable time, fol- 
lowed by the complete rupture occurring suddenly Note 
the sharp contrast with the results based on the conven- 
tional discrete-time recursive dynamics, where the strong 
stress brings about immediate failure [lfj. 

We also examine the effects of noise (T ^ 0), which 
disallows the nontrivial solution in the stationary state 
and leads to f c = 0. Nevertheless Fig. ^displays that the 
rupture time can be very long unless / is not large. In a 
real system T is nonzero but usually very small, and the 
complete rupture practically does not occur in the time 
scale of interest. For large /, the rupture is accelerated 
in the presence of noise, as shown in Fig. 21 

In summary, we have introduced realistic continuous- 
time dynamics for fiber bundles and investigated the be- 
havior of the system under stress. The dynamics has been 
generalized to include uncertainty due to impurities and 
environmental influences. In its presence the system has 
been found always to break eventually, reflecting the ir- 
reversible nature of breaking, whereas in its absence all 
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